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1. INTRODUCTION 

Attempts to use widely superconducting large scale systems give an impetus to 
comprehensive investigations on methods of cryogenic stabilization. Comparing the 
alternative procedures, it becomes clear, particularly for superconducting power 
transmission lines, that the advantageous approach consists in cooling supercon- 
ductors by forced circulation of supercritical helium /!/. There are several 
papers /2-7/ presenting experimental data on forced convection heat transfer to 
supercritical helium. It is worth while noting that these papers have successfully 
correlated their experimental data with an accuracy of about ±30%, but have failed 
to correlate the data presented in different papers. This is an evidence of either 
insufficient measurement accuracy or inadequate knowledge of the physical processes 
under study. Most experiments have been carried ouo with ducts of relatively small 
diameters (1-2 mm). As a consequence, the effect of geometry and body forces on 
heat transfer was negligible and, therefore, unaccounted for in the correlations. 

Ref. /8/ reported on the heat transfer to supercritical helium in vertical ducts of 
a comparatively large diameter (l8 mm); but for proposed designs of superconducting 
power transmission lines /9/, generally rather large (» 1*2 mm) horizontal cool- 
ant flow channels are envisioned. In the near-critical region, thermophysical 
properties are very sensitive to heat input; thus one would expect substantial sec- 
ondary flow effects in a heated long horizontal channel of sufficiently large diam- 
eter (much greater than lt2 mm). This assumption follows from the similarity 
theory; indeed, the considerable nonuniformity of wall temperature distribution 
over the periphery of a horizontal tube cooled by the flow of supercritical water or 
CO 2 has been verified experimentally and studied in depth /10-21/. But a dearth of 
experimental data on combined free and forced convective heat transfer to super- 
critical helium exists in literature. It is obvious that such information is of 
importance for superconducting power transmission line design when fault conditions 
are under consideration and in all cases when heat input to the superconductor becomes 
more than its steady state value - e.g. , due to leakage in vacuum space or short 
circuiting. The data may be also useful for designers of any large helium cooling 
systems. 

Along with the experimental approach there are attempts to obtain data on heat 
transfer to supercritical helium using the similarity theory in conjunction with 
generalized empirical relationships for other supercritical media /22/. The at- 
tempts appear to be rather successful /23/ but reliable results for an unstudied 
range of parameters with large thermogravitational effects could not be insured. 

In this paper we have extended our preliminary results /2k/ on the measurements 
of combined free and forced convective heat transfer to supercritical heliian flow 
and attempt to explain and correlate the data. 

The experimental conditions covered correspond to usually recommended pareun- 
eters for SPTL projects /9/. 

2. APPARATUS & ANALYSIS 

The experimental set-up is that described in /2k/ and is in part repeated here 
only as a convenience to the reader. 
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In the flow loop schematic. Fig. 1, helium flows from the refrigerator- 
liquifier 1 to the test channel 5 via a cryogenic pipeline 2 and distribution head 3 
of control unit J*. For closed loop operation, the helium flow from channel 5 is 
delivered through a return pipe to refrigerator 1 cooling the radiation shields of 
channel 5, unit and pipeline 2. For open loop operation, the helium flows from 
channel 5 to a gasholder via an open valve of unit U, heat exchanger 6 and flowrate 
measuring section 7. In both cases the main parameters of the flow (p, T, G) were 
determined by the refrigerator-liquifier setting and controlled by means of throttle 
valves of head 3 and unit 

Liquid nitrogen cooled shields were used to reduce heat input to the helium 
flow loop from room temperature. 

Fig. 2 shows the stainless steel horizontal circular test channel. A manganin 
wire heater is wound with a 4 mm pitch on the 2.85 m long, 2.0 cm o.d. , 1.9 cm i.d. 
heated section of channel 5* A stainless steel, 6 mm o.d., 5 ni® i.d. rotatable co- 
axial probe enables one to obtain temperature profiles at two cross-sections of the 
flow. When the probe is rotated, its thermocouples describe circles of ^4, 5.5, 7, 

8. 5 mm radii. 

The inlet and outlet 'I'^2 bulk fluid temperatures were measured with 

arsenide-gallium resistance thermometers introduced into the flow downstream of 
mixing chambers at either end of the test section where velocity and temperature 
profiles over the cross-section are averaged. Temperature profiles along the test 
section wall, over its periphery and at two flow cross-sections were measured by 
copper-iron thermocouples /25/. Positions of the thermocouples and thermometers are 
shown in Fig. 3. 

Without heating, and T^2 agreed to within 0, 02 K. Errors in temper- 

ature measured with differential thermocouples were less than 0.2 K. The variation 
of bulk fluid temperature T^ along the test section was calculated from inlet 
te.mperature , pressure and heat input (the pressure drop along the heated section was 
negligible). For open cycle operation the mass flowrate obtained from the heat- 
balance was compared with the more accurate orifice flowmeter (at room temperature); 
the comparison served as one validity criteria for the experimental procedure. 

For closed loop operation, the flowrate was determined only from the heat- 
balance equation but preliminary calibration with the orifice flowmeter was carried 
out for every run. Values of mass flowrate determined by these two methods agreed 
to within ±5^. 

Helium properties needed in the analysis were determined from /26/, stainless 
steel conductivity data / 27 / were used to calculate the inner wall temperature of 
the heated section by means of an iteration procedure. 


r 


Experimental results 

A more detailed analysis of the data than appears in /2li/ is presented as 
Table I; also listed are computed parameters required in the analysis that follows. 

The use of the equivalent-hydraulic diameter Djj as a characteristic size 
agrees well (±3?) with the correlation 


Nu 0.45 

Nu^ = - 2.1 + Pr \D / 

proposed / 29 / for annular tubes with heat input only to outer tube, for the range of 
parameters covered in these experiments (d/D = 0.3). 

Wall temperature distribution over the heated section periphery 

The temperature distribution (Fig. h) at x/L = 0.57 was measured by eight 
thermocouples equally spaced on the tube periphery. The profiles are typical of 
flows in heated horizontal tubes, that is, the highest temperatures are observed in 
the upper surface while the lowest ones occur on the lower surface. We compared 
profiles obtained (60 < Grq/GrL < 600; x/L = 0.57) with those calculated from cor- 
relation /31/ 




= cos 


( 1 ) 
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where for the steady fK'W with near-ideal properties at Re = 8.5 * 10^ to 5.2 * 10^ 
Grq = (0.2 to ..) • u''^L = 1.0*,' t-o 1.3 where Orq/Grp, < 250, n * f(Grq/GrL) and 

n = 2 at Orq/Gr^ > 250 to iOO (solid curves in Fig. u). It can be seen that ex- 
perimental and calculateil vahies a^iree qualitatively. But it sliould be noted that 
In all regimes (test im:is ) tt-.e measured temperature exceeds tiie calculated one at 
i = Ti/i* and in regimes 2, 6 and 12, the first temperature being even higher than 
that on the upper sia-face. It is unlikely that this could be attributed to experi- 
mental errors since each point was covered by two thermocouples. Further wall tem- 
perature ripples can be assuuied to be associated with features of flow structure in 
horizontal t\»bes aiid catised I’y Uie effect of secondary flows. 

Temperature distribution in iieliura flow 

Fig. 5 gives the typical temperature distribution obtained for regime 2 at 
x/L = 0,1. Isotheras are plotted on the basis of local flow temperature profiles 
measured by the probe. The pattern of Isotherms is typical of flows with mixed con- 
vection. It is of if.terest to note that the maximum difference between local tem- 
peratures in the upper and lower parts of the flow were measured to be o.l K while 
the maximum difference between measured wall temperature was about 3.0 K. This 
significant discrepancy can be attributed to the presence of low temperature zones 
in the flow. The isotherm patten: found evidence also of the existence of a cell- 
wise flow structure I rought al>out by secondary flows of thermogravitational origin. 
Examining the changes in the pattern of isotherms, it can be concluded that the flow 
includes six vortices. 

In the upper part uf t*.e channel the gradient dT/dh on all lines with h = 
const is directed towirl- the vertical axis. Thus, the isotherms indicate that 
there exists an up'ward fli w along the vertical axis. At the angle of * ti/ 3 the 
isotherms undergo change.*:; curves begin to bend downwards approaching the wall. 
Evidently, the flow heated up from the wall in the middle portion of the channel and 
moving along the wall downwaris and, as a re.sult of their interaction, the first 
flow begins to move towards the channel axis. Also, tliere is a zone with the rapid 
rise in temperature gradient, dT/d^ at $ * 0.8 it which indicates the existence 
of secondary flow. Julging by isotherm shapes, it may be assumed that the flow in 
the zone is direct.'d from the channel axis to its wall. L'ue to a relatively great 
distance between tliermocoui le Junctioi:s (about 1.5 m.Ti) and to low temperature gra- 
dients at 0 * ir, it prove! to I'e impossible to obtain the complete picture of sec- 
ondary flows at tlie bottom of t!ie cliaimel. However, the character of flows in ad- 
jacent regions of tlie cliunnel suggests tliat an upward flow takes place at iji * n 
along the vertical axis. ilierefore t !ie temperature distribution in tile flow imp^lies 
that there are at least six large vortices in the channel which are in contrast to 
the isothermal surfaces, constant property pliotographs in /2i*/. 

The flow pattern olitainel allows some features of temperature distribution over 
the wall periphery to be explained. Evidently, tlie iiighest wall temperature is un- 
likely to be found in tlie u[ per surface exposed to secondary flows coming from 
cooler fluid layers but near the region with = tt/3 where two secondary flows 

with the highest energies .merge. 

This is respons ible foi- a relative rise in wall temperature at = v/k. At the 
same time, it might be exjiected that the wall temperature at the lowest surface ex- 
ceed those at ^ = 0.8 n. But this appeared to be true only in regime 2 while in 
regime 15 the opposite t ffeet was observed and in all other regi.mes the feature was 
not noted at all. Thus, it cun be ussuiaed that in the lower portion the flow struc- 
ture is more complex than that shown in Fig. 5 and depends on regime parameters; per 
haps the behavior is more closely related to developing forced convection flow. 

Wall temperature distribution along tl:e heated section 

Fig. 6 shows wall temperature profiles along the upper and lower surfaces of 
the channel. Calculated values of bulk fluid temperature are also plotted. 

The naximiun temperature lUffr-rence over tlie wall periphery is seen to he equal to 
8.7 K (regime 10) whicii should be considered as a rather significant value if it is 
taken into account tliat bulk fluid temperature was 12.5 K. Also the data presented 
show that entrance and exit portions of the duct have some influence on the wall 
temperature profiles. Comparing changes in wall and flow temperatures it may be 
assumed that the lengths of the flow thermal stabilization are different for upper 
and lower wall surface that Is, the length for upper surface appears to be 1.5 times 
less than that for tlie lower one. 
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Heat transfer coefficients 


Assuming (l) adequately describes the peripherial surface temperature, and tak- 
ing n = 2 


^w 




where ^ S" ^’i('t')d4i; and s 




since Ku = — it follows that* Nu. = 

* t 




Nu for three tube cross-sections (x/L = 0.375; 0.57; 0.75Q) where 

iJi — TT 


thermal stabilisation of the flow has been achieved. 

An attempt was made to generalize the data by the relationship /23/ proposed for 
supercritical turbulent forced flows at large Reynolds numbers 




H f. 


Nu = 
r 


Nu 

Nu^^ 


f7==^=r-T H r, 

yy 0. a ii' + o.2 + iy 


(2a) 


(2b) 


where for our results the local parameter v*> replaces and fj^ or f 2 group 
the data with the same rapport, f^ f;; > f. As for the lower surface, our data 
agree with Eq. (2) satisfactorily enough; 73)* of points deviate from the predicted 
curve (Fig. 7(a)) only by i20‘*. 

Therefore, we believe tJie heat transfer on the lower wall surface within the 
range of parameters covered to be weakly dependent on the thermogravitational effect 
and, thus, can be estimated by usual correlations (Eq. (2)). 

The application of a single parameter to correlate experimental data ob- 
tained for the upper wall surface proved to be inadequate; all experimental points 
fall below the predicted curve by 5C‘* on the average. As a consequence, the corre- 
lation of Nu^ values by Eq. (2) was also unsuccessful, and experimental data were 
by 20 ! 30,* below the p'redicted curve, see also /2^/. This is quite conceivable 
since Eq. (2) does not include a term accounting for the effects of secondary flows 
due to thermogravitation on heat transfer. 

To correlate experimental data on the upper surface we modified the relation- 
ship /31/ 


Nu 




1 + 0.03'. 






•ftj 


0.0h5 - ^3 


where Nuip is the turbulent Nusselt number to: 

Nu 


Nu 


to 


= f • f.. 


( 3 ) 


Form similar to van der Waals binary mixture rule and differs from that of 
/2U/ where Nu^ = xNuq + (l - x)Nu^ and ^ = x^iq + (l - x)t(;^ with x » 


rm 


dt and D„ = 
n 


_i 

IT 


D-d. 
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This data representation as compared to (2) turned out to be more attractive, since 
nearly all points were within the range of ±30!5. Heat transfer deterior^ion at 
= 0 due to free convection influence causes the genera^reduction of Nu^ in com- 
parison with purely forced flow. The deterioration of Hu* seems to be related to 
the trend to suppress the turbulent flow in the upper portion of the channel. It is 
evident that the suppression exceeds additional flow turbulence near the lower 
surface. Using correlations (2) and (3) in conjunction with (l), the expression to 
estimate average heat transfer for combined free and forced convection for a hori- 
zontal channel can be derived 



where ti' = ip* = 1 6(~w^ -T^),Twa “ wall temperature averaged over its peri- 

phery', Eq. (1) has been substantiated by our experiments and predicts experimental 
data with an accuracy of about ±20' (Fig. 7). 

3. CONCLUSIONS 

An experimental study on heat transfer to low temperature helium flow in a 
horizontal channel where thermogravitational forces affect the flow considerably has 
been carried out. 

It has been found that there are six vortices causing an upward flow in the 
upper portion of the channel which gives rise to a shift in the location of the high- 
est temperature from the upper:; .st point on the wall surface to the point at 
* ir/3. 

Wall temperature profiles have been obtained. The contradiction between the 
measured temperature distribution over the tube wall periphery and that given by 
Eq. (l)) can be attributed to the effect of secondary flows caused by thenrogravita- 
tional forces. 

The thermal stabilization length on the lower wall is demonstrated to be on the 
average 1.5 times greater than that on the upper surface. 

Experimental data on local heat transfer both on the lower and upper wall sur- 
faces have been correlated using the parameter Grq/GrL, accounting for th^e effect 
of thermogravitational forces, and the local parameter = 1 + 6(Tw,j| - T^). It is 
found that the heat transfer coefficient for the lower surface has the same value as 
in the case of flow without free convection effects. 

To estimate the heat transfer coefficient averaged over the wall periphery, it 
is proposed to use for horizontal tubes the dependence Nur|^ = f(ijr^) where = 



1*. NOMENCLATURE 


X 

L 

D 


axial distance from 
heated section in- 
let 

Q 

nDL = Fp 

heater electric power 
heated surface area 

length of the b.eated 
section 

0 = _ 

3T 

isobaric expansion co- 
efficient 

(heated section); 
inner diameter 


gravitational acceler- 
ation 


d 



P 


(probe); outer 
diameter 

equivalent hydraulic 
diameter 

temperature 

viscosity 

heat conductivity 

mass flow rate 

specific heat 

pressure 


P 



density 

angle; = 0 uppermost 
tube surface; ^ * n 
lowest tube surface 

Reynolds number 
Prandtl number 
Grashof number 


I 

\ 
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' “ Nu 

(T - T, )X 
w^ b 

local Nusselt 
number 

Gr 

°"l 

thermogravitational forces parameter 
(see /13/, /30/, /31/) 

41 


h 

distance from the horizontal surface 

1 6(T^ - = 4- 

non- isothermal 
parameter /23/ 


F 

1 + 6(T - T ) = H; 

local non- 


P 

W D 9 

4* 

isothermal 




parameter 

Subscripts: 

0.023 Re°‘®Pr°'^=Nu^ 

Dittus-Boelter 

4 

local 

dependence 

C 

upper surface 

a 

heat transfer 

IT 

lower surface 

Nu 

- — ■ Nu 
% 

coefficient 

W 

wall 

reduced Nusselt 
number 

b 

bulk 

="l 

limit value of 
the Grashof 

- 

value averaged over the wall peri- 
phery 


number; 




Gr^_ = 3*10“^-Re^ 

•’5 [1.2 

(Pr2'3.„jp^0.5 
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L' ETUDE EXPERIMENTALE DE LA TRANSMISSION DE CHALEUR PAR CONVECTION AU HELIUM AUPRES 
DE CRITIQUE DANS LE CANAL HORIZONTAL 

SOMMAIRE: On mesure les coefficients de la transmission du chaleur au helium auprfes 

de critique, coulant dans un tube chaufe de longueur horizontale 2.85 “ et de 
diamfetre int^rieur 19 mm. La gamme des paramfetres exp6rimentaux comprend 6.5 - 15 K; 
1.2 - 3.0 bar aux flux de chaleur Jusqu'i 1000 W/m^ et les nombres de Reynolds sont 
de I'ordre de 9*103 h 2.105. On trouve une uniformity considerable dans la distribu- 
tion du coefficient de la transmission de chaleur sur la peripherie du tube. La dif- 
ference des temperatures entre les points sup€rieurs et inf6rieurs de la section 
transversale du tube, fait en I'acier inoxydable, etait figale h 9 K. On remarque 
que la temperature maximale de paroi sur la peripherie deplace relativement du point 
superieur du tube. Les mesures locales de la temperature dans la section trans- 
versale de flux ont montre la presence de I'ecoulement tourbillonnaire au canal. 

Le systfeme tuorbillonnaire obtenu dans le flux explique le deplacement du point de 
la temperature maximale de la paroi de tube relativement du point superieur. On 
propose la dependance pour les coefficients locals et moyens de la transmission du 
chaleur dans le tube horizontale. 
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Experimental results; heat transfer to fluid helium flowing 
through a horizontal tube 
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the test seotion. 


















